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Abstract 

A simple model for electromagnetic wave propagation through zero- 
temperature plasma is analyzed. Many of the complexities of the plasma 
state are present even under these idealized conditions, and a number of 
mathematical difficulties emerge. In particular, boundary value problems 
formulated on the basis of conventional electromagnetic theory turn out 
to be ill-posed in this context. However, conditions may be prescribed 
under which solutions to the Dirichlet problem exist in an appropriately 
weak sense. In addition to its physical interest, analysis of the cold plasma 
model illuminates generic difficulties in formulating and solving boundary 
value problems for mixed elliptic-hyperbolic partial differential equations. 



1 Introduction 

Among the many equations of mathematical physics vifhich change from el- 
liptic to hyperbolic type along a smooth curve, only the equations for tran- 
sonic flow have received sustained attention. In this brief review we consider 
elliptic-hyperbolic equations originating in a simple model for the propagation 
of electromagnetic waves through zero-temperature plasma. Solutions to such 
equations are likely to have significantly weaker regularity than solutions to the 
linearized equations of transonic flow. Recognizing the interdisciplinary nature 
of the topic, we assume a familiarity with physics but not necessarily plasma 
physics, and analysis but not necessarily elliptic-hyperbolic equations. However, 
the physics is confined to a review of fundamental results in Sec. [U whereas the 
mathematical results of Sec. [3] are somewhat more technical. There we consider 
the extent to which problems formulated primarily for linearized equations of 
gas dynamics possess analogies for equations arising from a different physical 
problem. Continuing such investigations in various physical and geometric con- 
texts (c.f. [29j). one may hope to obtain eventually a natural theory for linear 
elliptic-hyperbolic partial differential equations. 
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1.1 Physical background 

The plasma state is characterized by the dorainance of long-range, nonlinear 
effects. For matter in such a state, it is particularly difficult to obtain mathe- 
matical problems which can be stated with a satisfactory degree of rigor, and 
for which solutions can be shown to exist. Without a proof of the existence and 
uniqueness of solutions — which, in particular, specifies the function spaces in 
which solutions lie — it is hard to place appropriate boundary conditions on 
numerical experiments and to gauge the reliability of the results obtained. 

If one hopes to obtain a tractable mathematical problem, it is usually nec- 
essary to impose harsh assumptions on both the plasma and the applied field. 
Perhaps the harshest of these fixes the temperature of the plasma to be zero. 
This permits one to neglect altogether the fluid properties of the medium, which 
is then treated as a linear dielectric. Somewhat surprisingly, the assumption 
of zero plasma temperature is a useful first approximation to the products of 
tokamaks: low-density plasmas which are remarkably free of expected high- 
temperature phenomena such as collisions and wall effects. See the remarks 
in the introduction to [36] and the more detailed discussions in [39]. More 
generally, the cold plasma model approximates the effects of small-amplitude 
electromagnetic waves, propagating with phase velocities which are sufficiently 
large in comparison to the thermal velocity of the particles. 

We note that the term cold plasma is highly ambiguous. Although we take 
this to imply zero temperature, in the astrophysics literature interstellar plasmas 
on the order of 10'' K to 10^ K are typically referred to as "cold" (see, e.g., 
[TT|). Very recently, "ultracold" neutral plasmas, having electron temperatures 
in the range from 1 K to 10^ K and ion temperatures ranging from 10~'^ K 
to 1 have been created experimentally. The cold plasma model explored 
in this paper is apparently too simple to yield quantitative insight into those 
plasmas. In particular, the fluid dynamics aspects of experimental ultracold 
plasmas appear to be non- negligible (c./. Sec. 3 of [T7]). 

The other physical hypotheses imposed in this review are also quite restric- 
tive: Although the plasma is not assumed to be homogeneous, the inhomogene- 
ity is taken to be two-dimensional, so the governing equations for the model 
are also essentially two-dimensional. Moreover, the applied magnetic field is as- 
sumed in Sec. 12.41 to be longitudinal and the resulting wave motion confined to 
electrostatic oscillations. In Sec. I2.5l we consider electromagnetic waves, but we 
find (after reviewing a detailed analysis by Weitzner [39]) that elliptic- hyperbolic 
equations arising in the electrostatic case retain their validity as a qualitative 
model for the general case. 

For the most part, the outstanding mathematical problems relevant to the 
cold plasma model are boundary value problems for Maxwell's equations. The 
dielectric tensor for these equations will render them of elliptic type on one part 
of their domain and of hyperbolic type on the remainder, except for a smooth 
curve (the parabolic line) separating the two regions. Little is known about the 
formulation of well-posed boundary value problems for equations which change 
type in this way, especially as the equations that arise in the cold plasma model 
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appear to have certain fundamental differences from those that arise in gas 
dynamics. 

Careful reasoning about the mathematical properties of plasma models is 
not needed merely in order to prevent "mathematicians' nightmares." An ex- 
ample is known |26| in which the boundary conditions suggested by physical 
reasoning about the plasma lead to a mathematically ill-posed problem in the 
expected function space. Moreover, numerical experiments tend to confirm the 
difficulties that arise when the model equations are subjected to classical ana- 
lytic techniques; see |26| and various remarks in [39 . 

High-frequency waves can be modelled via geometrical optics. (Any propa- 
gating electromagnetic field will tend to have high frequency relative to the char- 
acteristic plasma frequencies; see, e.g., Sec. 2.4 of [23].) Mathematical problems 
that arise in the geometrical optics approximation are quite different from those 
that arise from applying Maxwell's equations directly, and we do not pursue the 
geometrical optics approach in this review. The complexity of the geometrical 
optics approximation is due to significant difference in magnitude among the 
terms of the plasma conductivity tensor at lower hybrid frequencies; see [34j 
and the references therein. 

The physics presented in Secs. l2.f Iand l2.2l essentiallv goes back to the work of 
Tonks and Langmuir [35] in the late 1920s. The results of Sec. 12.31 were already 
well known in the 1950s [1H[2S]; those of Sees. [23] and [23] date from the 1970s 
[2n][S3] and 1980s, respectively. In particular, Sec. l2.5l derives some fundamental 
analytic formulas introduced by Weitzner in |39j and [40] ; see also [16j . Section[3] 
is based on recent results, [271 [28] which extend analogous research on equations 
of Tricomi type — particularly |2T and [22] ; see also [41] , an earlier paper which 
is based on [24j . 

In the sequel, a subscripted variable denotes (usually partial) differentiation 
in the direction of the variable, whereas subscripted numbers denote components 
of a matrix, vector, or tensor. Differentiation of vector or matrix components 
in the direction of a variable is indicated by preceding the subscripted variable 
by a comma. Unless otherwise stated, a cartesian coordinate system is assumed 
in which the subscript 1 denotes a component projected onto the x-axis, the 
subscript 2 denotes a component projected onto the y-axis, and the subscript 3 
denotes a component projected onto the z-axis. In particular, x = {xi, X2, X3) = 

{x, y, z) and we denote the canonical cartesian basis by I z, j, fc ] . 



2 The cold plasma model 

A plasma is a fluid composed of electrons and one or more species of ions. 
Because it is a fluid, its evolution must satisfy the equations of fluid dynamics. 
But because the particles of the fluid are charged, they act as sources of an 
electromagnetic field, which is governed by Maxwell's equations. The presence 
of this intrinsic field leads to highly nonlinear behavior. Indeed, the dominance 
of long-range electromagnetic interactions over the short-range interatomic or 
intermolecular forces is often cited as the defining characteristic of the plasma 
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state. 

If the plasma is at zero temperature, then Amontons' Law imphes that the 
pressure term in the equations for fluid motion will also be zero, and the laws 
of fluid dynamics will enter only through the conservation laws for mass and 
momentum. In fact, because collisions can be neglected, the fluid aspect of the 
medium can be virtually ignored. The plasma is then represented as a static 
dielectric through which electromagnetic waves propagate. 

In particular, zero-order quantities — the plasma density, proportions of ions 
to electrons, and the background magnetic field — can all be considered static 
in time and uniform in space. First-order quantities — the electric field E and 
particle velocities v, are assumed to be expressible as plane waves: sinusoidal 
waves proportional to functions having the form exp [i (k • r — utt)] , where k is 
the propagation vector of the wave (not to be confused with the notation for 
the cartesian basis vector k); r is the radial coordinate in space; uj is angular 
frequency; = —1. Thus in cartesian coordinates, k • r = kix + + k^z. 

In the following we review basic elements of the physical theory that results 
from these assumptions. The material in Sees. I2.1H2.3I is completely standard 
and can be found in many sources. The classical reference is Ch. 1 of [33]; see 
also [T], [S], [13], and Sec. 2 of [ID]. More recent surveys include nodes 43-45 
of [12] and Ch. 2 of [37]. A recent review of theoretical plasma physics can be 
found in [6]. We employ SI units except where other units are specified. 



2.1 Equations of motion 

Consider a single particle of mass m, having charge q = ZSe, where Z is a 
positive integer, 6 equals 1 or —1, and e is the charge on an electron. Let the 
particle be subjected only to the Lorentz force 

Fi=g(E + vxB), 

where 

B = Bok. (2.1) 

Equation ()2.1|) implies that the applied magnetic field is longitudinal: its only 
nonzero component is directed along the positive z-axis. (In fact, there is little 
harm in assuming, somewhat more generally, that 

B = Bok + B {x, y, z) exp [i (k ■ r — ujt)] , 

with |B| « \Bo\.) 

The equation of motion for the particle is given by Newton's Second Law of 
Motion, that is, 

m^=F,. (2.2) 
In accordance with our assumption about first-order quantities, we write 
V {x, y, z,t) {x, y, z) exp [i (k ■ r - ujt)] , 
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or 

dv 

—— — —ILOV. 

dt 

Substituting this result into (|2.2p yields 

- imcov = q(E + v xB^ , (2.3) 

where 

E (x, y, z,t) =1] {x, y, z) exp [i (k • r - ut)] . (2.4) 
Initially we will take E to be a constant vector: 

E{x,y,z)^Eii + E2j + E3k (2.5) 

where Ei, E2, and E3 are constants, and similarly for v. 
Defining the cyclotron frequency 



n = 



qB, 



m 

Eq. (|2.3p has solutions v — {vi,V2,V3) satisfying 

^1 = , r n2^ (^^1 + ^^"^2) ; (2.6) 

TO (W^ — Si^) 

^2 = . (^^2 - i^r!£;i) ; (2.7) 

TO (OJ^ — iZ^j 

= -^E^. (2.8) 

TOW 

2.2 The dielectric tensor 

Although the above relations were derived for an individual particle, they also 
hold, in our simplified linear model, for each species of particle in a plasma 
consisting of electrons and iV — 1 species of ions. In particular, the plasma 
current can be written as the sum 

N 

] = ^n^q^v^, (2.9) 
i/=i 

where is the density of particles having charge magnitude \q,y\ = Z^e. 

In the sequel we will only consider the aggregate of particles, in which Eqs. 
(|2.ip - p.8p pertain with the quantities v, to, q, Z, S, and Q indexed by i^, where 
u = 1^ N. Introduce the electric displacement vector 

1 

D = vacuum displacement + plasma current ~ eqE H j, (2-10) 

w 

where eo is the permittivity of free space. It will be convenient to express (|2.10p 
in the form 

D = eoKE, (2.11) 
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where 

3 

D, = e^Y.^-J^3 (2.12) 

i=i 

and K — {Kij) is the dielectric tensor (also caUed the cold plasma conductivity 
tensor). The tensorial nature of this quantity reflects the anisotropy of the 
plasma due to the presence of an applied magnetic field. (Note that in the 
sequel the reader will be expected to distinguish between the notation K for the 
dielectric tensor, the notation Kij for the scalar element of its i*'* row and j^^ 
column, and the notation K, for the type-change function of an elliptic-hyperbolic 
equation.) 



Equations ((^ - (P?TT|) imply that 




(2.13) 



where s, d, and p are defined in terms of: 

i) the plasma frequency, which for particles of the v*-^ species is given by 



ii) the permittivities R or L of a right- or left-circularly polarized wave 
travelling in the direction k; these are given by 



N „2 



^2 Vo; + 5^D,i, 
u—1 ^ 



and 

In terms of these quantities. 



^ tt2 / 



and 



s = i(i? + L), 

d^\{R~V), 

v=l 

The mass of an electron is considerably smaller than the mass of any ion; so the 
squared ion cyclotron frequencies obtained from combining fractions in R and 
L can be neglected, leading to the approximate formulas 

N-l ^2 

^ \ \ ^ IIS (2 14) 

I/— 1 
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and 

N-l „2 
~ ^ w2 - wile + ^e^i ' y ■ ) 

In these formulas, the subscripted e denotes the value of the relevant quantity 
for the electrons and the subscripted ii, denotes that value for the i^*'* species 
of ion. Note that, by the same reasoning, the ion plasma frequencies can be 
neglected in the definition of p. 

2.3 The plasma dispersion relation 

The field equations for the system described in Sec. l2.2l are Maxwell's equations, 

VxE = - — , (2.16) 

VxB = AioU + eo^j, (2.17) 

where fiQ is the permeability of free space. 

From the form of Eq. (|2.4[) . it is clear that whenever E and B are plane 
waves, Eqs. (|2.16l) and (|2.17p reduce to the simpler form 

kxE = a;B (2.18) 

and 

k X B = — i^oj — w/^oEoE. (2-19) 
We can rewrite Eq. (|2.19p to read 

k X B = — cj^o f — + eoE 

= -w^toD = -eo^owKE, (2.20) 

where we have used (|2.10p and (|2.11l) in deriving the last identity. 

Now using (|2.18p . (|2.20p . and the elementary identity ^oEq — c~^, where c 
is the speed of light in vacuo, we obtain 

k X (k X E) = k X (wB) = a; (k X B) 



= -w^^oAioKE = - ( - ) KE, 



implying that 

k X (k X E) + KE = 0. (2.21) 

Define the index of refraction vector 



n — — k. 
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With this construction, the scalar index of refraction acquires a direction: that 
of the wave propagation vector k. In terms of n, Eq. (12.211) reads 

nx (nx E) + KE = 0. (2.22) 

Conventionahy, k (and thus n) hes in the xz-p\a.ne. Denote by 6 the angle 
subtended by the vectors k and B. Then (|2.22[) can be written as the matrix 
equation 



—id 

2 

p — n'^ sin^ 



9 

s — n-^ 




This matrix equation has a nontrivial solution precisely when the determinant of 
the 3x3 matrix vanishes. The condition for the vanishing of that determinant, 
the cold plasma dispersion relation is, geometrically, the equation for the wave- 
normal surface: 

An'^ - Bn^ + C = 0, (2.23) 

where the coefficients satisfy 

A = ssin^e + pcos^e, (2.24) 
B={s^ - d^) sin^ e + ps (1 + cos^ e) , (2.25) 

and 

C = p{s^ - d^) . (2.26) 

Because the left-hand side of Eq. (|2.23l) is a quadratic polynomial in n^, we 
obtain from the quadratic formula the solutions 

, B±F 

for F satisfying F'^ ^ B'^ - AAC. Using (|2:24l) - ((2:26t to write 
F^ ^ (RL - psf sin'' 6 + Ap^d^ cos^ 6*, 

we obtain 



{sv? - RL) (n2 - p) ■ 

These equations yield criteria for cutoffs where rt = 0, or resonance, where 
n — )■ oo. 

Physically, cutoffs and resonances correspond to a change in the behavior 
of the wave from possible propagation to evanescence. Mathematically, we will 
identify certain resonances with a change in type of the governing field equation, 
from hyperbolic (implying wave propagation) to elliptic (implying evanescence). 
These transitions may be accompanied, under certain conditions, by reflection 
and/or absorption of the wave. 
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Sufficient conditions for cutoff are p = 0, R — 0, or L — — tfiat is, a 
sufficient condition is C = 0. A sufficient condition for resonance is yl = 
wliidi, given Eq. (|2.24p . can be written 

tan^e = --. (2.27) 
s 

Particular cases of interest are 6 = (propagation parallel to the magnetic 
field) and = Tr/2 (propagation perpendicular to the magnetic field). We will 
be particularly interested in the hybrid resonances at 9 = 7r/2, which occur at 
frequencies for which s = 0. 

2.4 Electrostatic waves 

The electric field is said to be electrostatic if it approximately satisfies 

E==:-V$, (2.28) 

where <I> is a scalar potential. Equation (|2.28l) is satisfied locally by all time- 
independent electric fields and in an ordinary dielectric, the converse is also true. 
However in cold plasma there also exist time-dependent solutions of (|2.28p . Cold 
plasma has been characterized as a linear dielectric through which electromag- 
netic waves propagate. Thus these waves include, in distinction to ordinary 
dielectrics, the special case of propagating electrostatic waves. 
We write $ in the form 

<^{x, y, z;t) = LP {x, y, z) exp [k • r - iujt] , 

and add to Maxwell's equations (I2.16p . (|2.17p the additional equation 

divB^O, (2.29) 

which follows from Gauss' law for electricity. 
Equation (|2.28p implies immediately that 

V X E = 0. (2.30) 

This is most easily seen if we use differential forms, and note that in terms of 
the exterior derivative, E = d^, so ()2.30p is just the well-known property that 

dE = = 0. 

(Here and below we will switch from vectors to forms whenever the calculation 
is made more transparent thereby; but we will not change notation for the 
underlying geometric object, making the convention that the argument of the 
exterior derivative is always assumed to be a differential form.) In either the 
vector or form notation, identity (I2.30p follows from the equality of mixed partial 
derivatives. Applying the arguments relating (|2.16p to (|2.18p and (|2.17p to (|2.19p 
(translation into Fourier modes), we rewrite (|2.30p in the form 

k X E = 0. 
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This implies, by the properties of the cross product, the geometric fact that the 
vectors k and E are parallel We say that electrostatic waves are longitudinal. 
Physically, they appear as oscillations along the axis of the magnetic field. 

Thus we conclude that transverse waves, which propagate in a direction 
perpendicular to the magnetic field, must satisfy 

kE = 0. 

Again, differential forms are illuminating: The above identity becomes Sd^ = 0, 
where 6 is the formal adjoint of the exterior derivative d. But $ is a 0-form, so 
(5$ = automatically, and we find that transverse waves satisfy 

M$ + d6^ = A$ 0, 

that is, transverse waves are necessarily harmonic. 

2.4.1 Plane-layered media 

If we allow a plane- layered inhomogeneous medium (parameterized by x), the 
electrostatic potential has the form 

* {x, y, z) = ip{x) exp [i {k2y + hz)] , 

where kj is the j*'' component of the wave vector k for j — 1, 2, 3. Substitution 
of this form for the electric potential into Eq. (j2.29p yields, using (|2.12p . the 
single scalar equation [20] 

Kiiipx^ + {Kii^a: + icTo) = 0, (2-31) 

where 

CTO = fcs iKi3 + K3i) + k2 {Ki2 +K21), 

and zero-order terms in ip have been neglected. This equation has a power-series 
solution except where Kn vanishes. 

Explicit solutions of the model equation (|2.3ip under various physical as- 
sumptions are given in Sec. 1 of [33 , the Appendix to [20], and Sec. C of [TB]. 

It is easy to believe that inhomogeneities may develop in a plasma. For 
example, if the temperature is not exactly zero, the difference in velocity be- 
tween electrons and ions can be expected to destabilize an initially homogeneous 
distribution. However, it is difficult to imagine a force that will restrict these 
inhomogeneities to a 1-parameter foliation, which would be necessary in order 
to arrive at Eq. (I2.3ip . Formally, an electromagnetic potential leading to Eq. 
(|2.31l) could be induced by applying a driving potential to the metallic plates of 
a condenser. But in practice, this plasma geometry has little application either 
in the laboratory or in nature. 
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2.4.2 A two-dimensional inhomogeneity 

Suppose instead that the medium is a cold, anisotropic plasma with a two- 
dimensional inhomogeneity parameterized by two variables, x and z. Then the 
field potential has the form 

$ (x, y,z) ^ ip {x, z) exp [ik2v] ■ 

The electric field E is then given by 

E = -V$ = {E^,E2,E^) = - {ip.,e'^-y,ik2^e^^\v.e''''y) . 

Maxwell's equations for the electric displacement vector D = {Di,D2,D^) 
take the form 

= V-D=I?i,,+i?2,y+i?3,.- (2.32) 

We continue to neglect those terms which do not contain derivatives of as 
is assumed to oscillate rapidly. 

Because neither ip nor Kij have any dependence on y, the problem is effec- 
tively two-dimensional. Applying Eq. (|2.12p . the surviving terms of Eq. (|2.32p 
are (setting eg equal to 1) 

Di^x = - [KiiLpxx + Kii^a;Vx + Ki2(pxik2 + Kisy^zx + Ki3^xfz] e'^^y-, 

D2^y = - [K2i(Pxik2 + K23'Pzik2] e'^^y; 

■03,2 = - [K^ilfixz + K^i^z^x + K32ik2(pz + Ks^ifizz + Kss.^z'Pz] e*''^^. 

Collecting terms, we find that [33 

Kllipxx + 2cr(y5a;z + K^alfizz + a^x + Ct2ipz = 0, (2.33) 

where 

ai = Kii^x + ik2 {Ki2 + K21) + i^ai.z; 
a2 = Ki3^x + ik2 {K2Z + ^^32) + K^^.z- 

Two-dimensional inhomogeneities of the kind represented by Eq. (j2.33p can 
be expected to arise in toroidal fields, such as those created in tokamaks. 

The entries of the matrix K under our assumptions on Bq imply that cr = 0, 
so we can write Eq. (|2.33p in the form 

KiWxx + K3j,ifzz + lower-order terms = 0. (2.34) 

Equation (|2.34p is of either elliptic or hyperbolic type, depending on whether 
the sign of the product 
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is, respectively, positive or negative. 

The sign of Kn changes at the cyclotron resonances — il^. The cold 
plasma model breaks down at these resonances, as three terms of the dielectric 
tensor become infinite. The sign of Kn also changes at the hybrid resonances, 
at which 

(These resonances, which have both a low-frequency and a high-frequency so- 
lution, are hybrid in that they involve both plasma and cyclotron frequencies.) 
In particular, the sign changes at the lower hybrid resonance, 

IP W 

1 + 4 = - (2.37) 

where as before, the subscript e denotes electron frequency, and the subscript 
i denotes ion frequency. At the hybrid resonance frequencies the cold plasma 
model retains its validity. 

The sign of if 33 changes on the surface 

1 = E^' (2.38) 

the resonance at which the frequency of the applied wave equals the plasma 
frequency of the medium. We may suppose that an electromagnetic wave prop- 
agating through a plasma does so at a much higher frequency than any of the 
characteristic frequencies of the plasma. Otherwise, the plasma magnetic field 
would prevent the waves from propagating very far (c./. |18|). Thus in evaluat- 
ing (|2.35p and in the sequel we will take if 33 to be strictly positive. 

Borrowing the terminology of fluid dynamics, we will refer to resonances 
such as (|2.36p ~ (|2.38p as sonic conditions on Eq. (|2.34p . 

2.4.3 The type of the governing equation 

In order to understand the possible variants of Eq. (I2.34p . we consider the 
coordinate transformation (a;, z) — )■ (^ (x, z) , ry {x, z)) , where 

C = Kii (x, z) . 

In these coordinates, the higher-order terms of Eq. (|2.34|) assume the form 

ifii<^,, + if33'P.. = (ee' + i^saCf) m + 
{iixTla^ + K^i^zVz) nv + i^vl + K^iVl) Vm- (2-39) 

In order that the transformation [x, z) — > (^, rf) be nonsingular, we require that 
its Jacobian be nonvanishing, i.e., 

LVz - LVx 7^ 0. (2.40) 
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Because we want the coefficients of the cross term (^^^ to be zero in the new 
coordinates, we impose the condition that 

^^ccrh + K33^zV,^0. (2.41) 

Both ^ and are given, and it is easy for the two first derivatives of 77 to 
satisfy (I2.40p and (|2.4ip simultaneously. 
Two possibilities arise. Either 

i) and never vanish simultaneously, or 

ii) there exist one or more points (x, z) on the domain at which 

e(x,z)=6(x,z)-0. (2.42) 

In case i), the condition f = implies, via (|2.4ip and the assumption that 
K33 is positive, the accompanying condition r]z = 0. But if ^ and r]z both vanish, 
then the coefficient of (p,,,, in ()2.39|) vanishes; that is, 

^Vl + = 

whenever ^ = 0. Again using (|2.41l) . we obtain from Eqs. (|2.34p and (|2.39p an 
equation with higher-order terms having the form 

The denominator in the coefficient of cannot be zero: ^ and cannot vanish 
simultaneously, and if vanishes, then must be nonzero in order to preserve 
condition (j^^ . So Eq. (1^^ is of the form 

(2-44) 

where /C (^, 77) = if and only if ^ = 0, an equation of Tricomi type. 

In case ii), condition p.40p prevents rjz from vanishing when ^z vanishes. 
Thus if <^ and £,z vanish together, the coefficient of ipriri in (|2.39p will not vanish 
at that point. Thus in case ii) we obtain from (|2.34p . (|2.39p . and (|2.41l) an 
equation with higher-order terms having the form 

^Vl+K^svf''^^''^'' (2.45) 

where the numerator in the coefficient of (p^^, but not the denominator, is zero 
whenever ^ is zero. That is, Eq. (|2.45p is an equation of the form 

ICi^,v)m+^vv^O^ (2-46) 

where IC (^, ry) = if and only if ^ = 0, an equation of Keldysh type. 

See Sec. 1.2 of [5], [9], Sec. 1 of [26], and Eqs. (75)-(78) of [39] for arguments 
of this kind. 
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The regularity of weak solutions to equations of Tricomi type can be estab- 
lished by microlocal arguments; see [M] and |T5] and, especially, [30] and [3T] . 
These arguments appear to fail for equations of Keldysh type, and one expects 
weaker regularity for weak solutions to such equations. 

The crucial question is: does condition (I2.42[) occur in our physical model? 
The answer to that question is "yes." 

2.4.4 Geometry of the resonance curve {after Piliya and Fedorov) 

Returning to our original xz-coordinates, we set the elements Kn and K22 of 
the dielectric tensor equal to /C, and the element K33 equal to rj. The coefficients 
of the only other nonzero elements, K12 and K21, are zero in Eq. (|2.34l) . so only 
Kii and play a direct role in the analysis. The sonic condition is equivalent 
to the alternative: 

/C = (2.47) 

or 

JCsm^e + ricos^e = 0, (2.48) 
where (x, z) is the angle between the direction of Bq and the xz-plane; c.f. 

The singular points on the sonic line (|2.47p are the points at which this 
curve (which is not a generally a line in standard coordinates) is tangent to the 
projection of the force lines of Bq in the xz-plane — that is, the flux lines of the 
magnetic field. The singular points of the graph F of Eq. (|2.48p are the points 
at which the fiux lines of Bq are normal to F. 

This motivates the placement of the origin at a singular point of the sonic 
line, with the z-axis (the axis along which Bq is directed) tangent to the sonic 
line. The x-axis is directed along the inward normal to the sonic line, relative to 
the hyperbolic region of Eq. (|2.33|) . Then Ku and a both vanish at the origin. 
Taking both x and z to be small, one can write 

T 7 

= - + - (2.49) 
a b 

and 

K33 = -r,o, (2.50) 
where 770 is a positive constant. Scale x and z, via 

a; — > i = x/a (2-51) 

and 

z z = z I a^/r^, (2.52) 

in order to obtain dimensionless variables i, and z. In this way, one obtains in 
place of (|2.33|) the equation 

- (i + Az^) (pjj + (^25 - (y3i = 0, (2.53) 

where A is a constant, for the simple case in which the coefficients of cross terms 
vanish identically |33] . 
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2.5 Analytic difficulties in the electromagnetic case {after 
H. Weitzner) 



In this section we suppose that the electric field satisfies Eqs. (|2.4p and (|2.5p . 
but no longer assume that the electric field satisfies condition (|2.28p . Closely fol- 
lowing [39] , we attempt to study the resulting field equations using conventional 
analytic tools, in order to see what difficulties arise. 

Repeating the calculations of Eqs. (|2.16p ~ (|2.2ip in greater detail, we compute 

V X (V X E) = V X ( --^ j = V X {iujB) = 

iuj (V X B) iuj /io M + eo-^ 
iuj^oj + ii^Mo^o (— iwE) = iujfJ-oi + w^/xoeoE 
..„.'f',..EU.o.'D.„...^KE. (.34, 



Now 



V X E = (Ea^y - E2,^) I + (E^.,, ~ E^^,) j + (Sa.x - Ei,y) k. 



It is obvious that this quantity vanishes identically in the electrostatic case: 
apply p.28p and equate mixed partial derivatives. But if V x E is itself a 
gradient, then the quantity V x (V x E) vanishes for the general case as well. 
We will understand the seriousness of this latter difficulty once we evaluate the 
left-hand side of Eq. (|234)) . Exphcitly, 

V X (V X E) = {E2,xy - Ei^yy - Ei^,, + E^^^^) i+ 

(Es^yZ ~ E2^ZZ — E2.XX + El y,j.) j + {El zX — E^ r^x — E'i yy -^2,21/) k. 

Applying (|2.4p . (|2.5p to the right-hand side, we obtain the algebraic expression 
[kik2E2 - (fci -I- kl) El + kik^E^] I + [k2k3E3 - {kj + kj) E2 + k2iEi] j 

+ [ksiEi - {kl + kl) E3 + A:32£;2] k. (2.55) 
This object can be written as the matrix operator 



iE : 



The system p.54p is uniquely solvable if and only if the operator L can be 
inverted — that is, if and only if 

- {kl + kl) kik2 fcifca 
det ( k2ki - {kj + kf) k2k3 | ^ 0. 
k^ki k3k2 - (kf + kl) 



{kl + /Cg) 


kik2 


kiks \ 


(E, 


fc2^1 


-{kl + kf) 




E2 


k^ki 


kzk2 


-{kf + kl) J 


V E3 
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But it is easy to check that this determinant vanishes identically for all {ki ,^2,^:3) . 
Of course (|2.55p is just a translation of V x (V x E) into Fourier mode. Because 
the symbol of a differential operator is a natural generalization of the idea of 
Fourier modes, we can interpret the foregoing computation to mean that the 
symbol of the differential operator V x (Vx ) vanishes identically. This is a 
serious obstacle to understanding (|2.54l) . As Weitzner notes in [30] , neither the 
type of Eq. (|2.54p (which is given by the sign of the symbol) nor the order of 
the equation (which is given by the degree of the symbol) are determined by 
standard analytic methods. 

2.5.1 Choices of potential and gauge 

It is therefore necessary to impose an additional hypothesis. A natural one is 
that an electromagnetic potential exists. But in distinction to the electrostatic 
case, we do not assume that E can be derived by simply taking the negative 
gradient of a scalar field. 

In order to compare our computations with the extensive expositions in 
and [10] we adopt, only in Sees. 12.5.11 and 12.5.21 the convention that the 
time- harmonic dependence is of the form exp [iujt] in units of c/uj. (This is in 
distinction to (j2.4p .') Because our time derivatives usually end up being taken 
twice, this only has an effect on the sign in a few intermediate calculations. 
However, with this sign convention. Maxwell's equations for plane waves assume 
the slightly different form 

V X E = iB, (2.56) 

V X B = iD = iKE. (2.57) 

The first choice of potentials is to let the vector A denote the magnetic 
potential and to introduce a second, scalar potential, $. We then write 

B = V X A (2.58) 

and 

E = -iA-V$. (2.59) 
Taking the curl of the second equation, we obtain 

V X E = -iV X A- V X V$. (2.60) 

Evaluating the last term on the right-hand side of (|2.60p using differential forms, 
$ is a zero-form and 

VxV$ = d2$==0. (2.61) 

Equations (|2.58p and (|2.6ip imply that ()2.56p is satisfied under condition (|2.59p 
for any smooth choice of A and $. 

Notice that we automatically obtain from hypothesis (|2.58p an extra condi- 
tion 

V • B = V • (V X A) , 
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which is to say, in terms of differential forms, that the 2-form B and the 1-form 
A satisfy 

dB = (fA ^ 0. 

In order to evahiate (j2.57p . we notice, continuing to interpret the magnetic 
potential A as a 1-form and $ as a zero-form, that if g is defined by 

5 (A) = A + df, 

where / is a smooth 0-form, then 

d {g (A)) = d{A + df) = dA + d^f = dA = B, 

so the magnetic field remains invariant under the gauge transformation g. More- 
over, because Sf = for any zero-form /, we have 

Af ^ ~{Sd + dS) f ^ ~6df. 

Thus, given any smooth potential A, we can choose / to satisfy the Poisson 
equation 

A/ - 5A, 

in which case 

S [g (A)) =5{A + df) = 5A- Af = 0. 

We call the gauge produced by such a g a Coulomb (transverse, radiation, or 
Hodge) gauge. In vector notation, 

iV ■ g (A) = 0. 

Computing (|2.57|) in the Coulomb gauge, we obtain [39] 

Ag (A) - iKV$ + Kg (A) = 0. 

Computing the symbol a of this operator by the same method that was applied 
to the operator V x (V x E) , we find that cr = — Ikj"* (Kk) • k, a polynomial of 
degree six in k. That the corresponding system is of order six is an expected 
result for a system of two first-order equations for vectors in M"^. 

Replacing the Coulomb gauge by a slightly more complicated gauge in which 

iV ■ (K*A) ^ 0, 

where the superscripted asterisk denotes the adjoint matrix, we obtain a self- 
adjoint operator in (j2.57p . This more complicated gauge can be constructed 
by the same general method that led to the Coulomb gauge, provided that we 
solve a slightly more complicated Poisson problem. The symbol corresponding 
to this self-adjoint operator can also be calculated by the methods introduced 
earlier, and that symbol is also a sixth-degree polynomial in k. 
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However, we can obtain a fourth-order system, which is more convenient for 
analysis, if we impose an additional hypothesis: that the plasma has axisymmet- 
ric geometry. While this is a very strong hypothesis, it is satisfied by plasmas 
produced in tokamaks. 

In order to motivate the choice of potential in this case, we make a few 
preliminary calculations. Only for the remainder of this section, the subscripts 
r, 9, and z when affixed to a vector are to be interpreted as the radial, angular, 
and axial components of the vector unless preceded by a comma; if preceded by 
a comma, they are to be considered partial derivatives in the direction of the 
component. (The subscripted- variable notation for partial derivatives of scalar 
functions remains unchanged.) Adopting the basis 

Uj. = cos 9i + sin 6j, 

ug — — sin 9i + cos 0j, 
Uz = -k, 

we recall that in the axisymmetric case we can write 

E = {Er (r, z) Ur + Eg (r, z) rug + Ez (r, z) u^) e 



and similarly for B. If m — 0, the waves preserve the axisymmetry of the 
underlying static plasma medium, as the wave vector satisfies k = {kr,0,kz) ■ 
We will restrict our attention to this simple special case, in which 

V X E = i [Ez,gUr + Er.z {rug) + {rEg)^ Uz] 
r 

-- [Ez.r (rug) + rEg,zUr + Er.gU^] = 

r ' 

-Eg^zUr + {Er,z " E^.r) Ug + - [rEg) ^ U,. (2.62) 

r 

Thus (|2.56p implies in particular that 

- iEg^z = Br (2.63) 

and 

{rEg),r 



-iB^. (2.64) 
r 

Just as Eqs. (|2.63p and (|2.64l) imply, using (j2.56p . that Br and B^ can each be 
expressed in terms of derivatives of Eg , so it is possible to use (|2.57l) to show that 
the other cylindrical components of E and B can be expressed as appropriate 
derivatives of Eg and Bg. This will allow the angular components of E and B 
to play the role of potentials for the two fields. 

Applying (|2.56l) to the middle term of the last identity in (|2.62p yields 

Er,. - E,^r = -iBg. (2.65) 
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Because E and B have exactly analogous forms and the left-hand and middle 
terms of (|2.57|) is exactly analogous to (|2.56p with a change of sign, we can 
immediately write 

and 




r 



Now the extreme right-hand side of (|2.57p yields Er and (see Eqs. (22), (23) 
of [32]) and one obtains 

Br,z - B,^r ^iDg^i {KgrEr + KggEg + Kg.E,) , (2.66) 

completing the system of equations for Eg and Bg. 



2.5.2 Variational interpretation 

Continuing to adopt the special hypotheses and special notation of Sec. I2.5.1[ 
we continue to review the analysis in |39| of geometry-preserving plane waves 
in an axisymmetric plasma. 

Equations (|2.65l) . (|2.66p can be associated to an energy functional: 

£ = J{[W {rE*g) ■ V {rEg)] /r^ + [V {rB*g) ■ KV {rBg)] /r^A 
[rBD ^ {K„Krg - K.gKrr) It + Bl , {K^gK,, ~ K,,kA /A 



-^iEg 

~lE*g 

where 
and 



{rBg) ^ [K.rKer ~ Kg.Krr) /r + Be,, (KerK, 



/A 



-B;Be ~ E*gEe [det (K)] /A}r drdz, (2.67) 



V 



A - KrrK, 



K K 



Provided K can be made self-adjoint, so can £. Form a right-handed orthogonal 
set (v, 0, u) , where 

u = sin j5r + cos /3fc 



and 



V = cos j3r — sin j3k. 



The basis is to be chosen so that u lies in the poloidal direction and v lies 
orthogonal to it; so we write the magnetic field vector in the form 

B = _Bo [cos aO + sin a (sin /3r -I- cos /3z)] , 

where a, /?, and Bq depend only on r and z. In this notation, the variational 
equations of £ form a second-order system in which the differential operator for 
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one of the equations is essentially the Laplacian £. We ignore that equation, as 
standard analytic methods can be applied to it. The differential operator for 
the other equation looks like 



(2.68) 



and that is the equation — in particular, the second of the two differential 
operators in that equation — that we will study in the remainder of this review. 
The term (|2.68p in Eq. (|2.67p can be written explicitly, in terms of the chosen 
basis, in the form 



rV 



V (rBg) 



rV 



C sin a 



V {rBg) X V 



{u-V){rBg)u 
fi cos a 



rA 



Bg 



(rA) ' [^i{C-Osinau-V{rEg) 
where 

c = e- 



i {f-i^ - ^C) sin a cos av • V (rEg)] , (2.69) 



N 

E 



and 



N 

E 



Equation (|2.69l) is only elliptic for negative values of Physically, this is the 
condition for so-called lower-hybrid frequencies, at which 



1 



c.f. (j2.37l) . Noticing that ^ is a function of r and z, define a new variable ry (r, z) 
so the curves ry = constant are orthogonal to the curves ^ = constant. Rewriting 
(|2.69p in (^, 77)-coordinates, the behavior of the solution depends on whether or 
not 

u • Ve = 0. 

This identity implies that flux surfaces coincide with resonance surfaces. In that 
case, Eq. (|2.69p is analogous to Eq. (|2.46p of Sec. 12.4.31 and is not of Tricomi 
type. The second-order terms of that equation can be written in the form 



L(u) = /(^,,7)Ku55 + Af(C,77)' 



(2.70) 



where u — u{^, rf) is a scalar function; / and M are given well behaved scalar 
functions near ^ = and, in addition, M is positive. 
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The physical model allows two further alternatives: If the curve representing 
the flux surface in two dimensions is collinear with the resonance curve as in 
(|2.70p . then the plasma can be treated as a perpendicularly stratified medium, 
which is essentially the case considered in Sec. 12.4.11 If the resonance curve 
is tangent to the curve representing the flux surface, then we are in a more 
mathematically and physically interesting case. In this latter case, the simplest 
model for the operator L of (|2.69|) is an operator for which the highest-order 
terms have the form 

L{u) ^ (x - y'^) + Uyy. (2.71) 

Note that this operator is closely related to the differential operator of Eq. (|2.53p . 
The two operators can be made virtually identical by replacing the coordinate 
transformation (|2.51l) by 

x^x = -x/a. (2.72) 



2.6 A conjecture about the singular set 

Methods for deriving the smoothness of solutions to the Tricomi equation appear 
to fail for an operator of the form ()2.7ip whenever the domain includes the origin 
of coordinates. This suggests the existence of a singular point at the origin, a 
conjecture which is supported by an analysis of characteristic lines. 

In order for a characteristic line to pass through the origin, the point(x, y) 
would need to satisfy the identity 

X = Ay2 (2.73) 

for some constant A, and also the characteristic equation for (|2.71l) . Substituting 
l|2.73p into the characteristic equation 

{x - y^) dy^ + dx^ = 0, (2.74) 

one obtains the equation 

^^3^^ ^ (1-A)y2' (2.75) 

or 

4A2 + A - 1 = 0. 

This polynomial has two real solutions; considering that the characteristic equa- 
tion (|2.75|) has two roots, one concludes [26] that four characteristic lines must 
pass through the origin — two more than pass through any other hyperbolic 
point. This motivates the suspicion that solutions of at least the equation 
Lu = will tend to be singular at the origin. It has been observed that an 
energy sink or plasma heating zone might be associated with such a singularity; 
see [T^ , [35] , [33] , [53] , and [10] for details on this and other issues raised in this 
section. 
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3 Analysis of the model equation 



Physical reasoning suggests that the closed Dirichlet problem, in which data are 
prescribed along the entire boundary of the domain, should be well-posed for 
the cold plasma model. However, the closed Dirichlet problem has been shown 
to be ill-posed, in the classical sense, for the equation 

(X - y"^) U^x + Uyy -f iuj; = 

on a typical domain [53]. This leads us to ask whether a well-posed problem 
with closed boundary data can be formulated in a suitably weak sense. In this 
section we address the "existence" part of that question. 

Because the operator introduced in Eq. ()2.71|) is not of Tricomi type at 
the origin, where it satisfies a condition of the form (|2.42|) . we expect weaker 
regularity than we have for operators which are uniformly of Tricomi type. In 
particular, although we can show the existence of very weak solutions in L^, 
we do not expect regularity for the closed Dirichlet problem. This lack of 
optimism is supported by numerical experiments |26| . Moreover, our methods 
are insufficient to determine the uniqueness of a solution. 

Denote by a bounded, connected domain of having piecewise smooth 
boundary 517, oriented in a counterclockwise direction; the domain includes 
both an arc of the sonic curve and the origin of coordinates in . (This insures 
that our equation will be elliptic-hyperbolic but not equivalent to an equation 
of Tricomi type.) 

Define |21j, for a given function /C {x,y) , the space (il; \JC\) and its 
dual. These spaces consist, respectively, of functions u for which the norm 

llullL^cn;!^!) = (^J \JC\u'^dxdy^ 
is finite, and functions u E L'^ {il) for which the norm 

l|w||L2(n;|K|-i) = (// m^^u'^dxdy^ 

is finite. Standard arguments allow us to define the space Hp (^5^) the 
closure of C§°{^1) with respect to the norm 



1/2 



\\u\\H^{n;K) = 
The i7g (ri; /C)-norm has the explicit form 



1/2 

[\IC\u^ + Uy + u^) dxdy 



1/2 

{\K.\ul +My) dxdy 



(3.1) 



(3.2) 



which can be derived from (jS.ll) via a weighted Poincare inequality. 

In the following we denote by C generic positive constants, the value of which 
may change from line to line. 
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3.1 The closed Dirichlet problem for distribution solutions 

Consider the equation 

Lu = /, (3.3) 
where / is a given, sufBciently smooth function of (a;, y) and 

^ ^ ^ dx^ ^ dy'^ ^ dx ^ ^ 

for a given constant k. By a distribution solution of equations (13.31) . p.4p with 
the boundary condition 

u{x,y) = Qy{x,y) dVL (3.5) 

we mean a function u e _L^(f7) such that G Hl{Vl] K.) for which L*^ e -^^(^^), 
we have 

(?.,L*0 = (/,0- (3.6) 

Here L* is the adjoint operator; ( , ) denotes the inner product on fi; ( , ) is 
the duality bracket associated to the norm 

II II 

iMH-Hn-X)^ sup — . 

Such solutions have also been called weak; c.f. Eq. (2.13) of [1], Sec. II. 2. In 
fact they are a little smoother than generic distribution solutions, as they lie in 
a classical function space. 

The existence of solutions to boundary value problems can be shown to 
follow from energy inequalities having the general form 

\\v\\v <C\\L*v\\u, 

where U and V are suitable function spaces. We will combine such an inequality 
with the Riesz Representation Theorem to prove the existence of distribution 
solutions; see for a general reference. 

Lemma 3.1. ( [28]). The inequality 

M\Hi(n;K) < C'll-^"llL2(n), 

is satisfied for u e Co(ri), where the positive constant C depends on il and K.; 
L is defined by \S.4^ with k G [0, 2] ; K. = x — y^. 

Proof. We outline the proof; for details, see [5S], Sec. 2. Initially, let 1 < k < 2, 
and let (5 be a small, positive constant. Define an operator M by the identity 

Mu ^ au + bux + cUy 
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for a = —1, c = 2 (2(5 — 1) y, and 



exp(2(5/C/Qi) ii{x,y)en+ 

exp(6(5/C/Q2) if(x,2/)er2- ' 

where 17+ = {(a;, y) <E | /C > 0} and fl^ — Choose Qi = exp (2(5/ii) , 

where /ii = max 



/C. Define the negative number 112 — min, s^pr^/C and 



let = exp (/X2) . Notice that & < Qi on f2+ and 6 > (^2 on fi" 

We will estimate the quantity {Mu, Lu) from above and below. As in the 
Tricomi case [Slj, one of the coefficients in Mu fails to be continuously differ- 
entiable on all of $7. When integrating this quantity, a cut should be introduced 
along the line /C = 0. The boundary integrals involving a, h, and c on either side 
of this line will cancel. 

The boundary terms vanish by the compact support of u. Integration by 
parts yields the identity 



[Mu, Lu) ^ 



au^ + 2f5uxUy + ^ui dxdy 



where 



for 



and 



si+uo- 

Cy bx 



K. 



b- cy, 



a\n- 



Qi 



6IC + 2{l-2S)y^ 



3— -2) 5|/C| +2(1-25)2/2 



/3=-[c(«-l) 



by] 



lib. 



y[2S{b/Q,) + {K-l){2S-l)] < \y\ 
y[6S{b/Q2) + (K-l){2S-l)]<n\y\ 

2{l-S) + S{b/Qi) mn+ 
2(1- S) + 36 (6/(32) in n- 



in Q+ 

in n- ' 



On r2+, for any scalars ^ and 77, we have 

aS,^ + 2/3^7] + 7772 > _ (-^2^2 ^ ^2^ _^ ^^2 ^ 



0^ 



(1 - 26) + 



6IC + {l~A6)y^ + [k- 

> 6 {jce + v') , 

provided 6 is sufficiently small. On 17^, 

a£,^ + 2(3^Tj + 7772 > a^f - 2 (y^^^ _^ ^2^ ^ ^^2 ^ 



Q2 



2 6\IC\ - 45y2 + K 



Q2 
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>5{\K\e + v^). 

Arguing in this way on each subdomain (and taking = u^, r] = Uy), we obtain 

{Mu,Lu)>S\\u\\l,^^.,^y (3.7) 

For the upper estimate, we have pTj 

iMu,Lu)<\\Mu\\L2\\Lu\\L2 < C iK,n)\\u\\Hif^^.ic)\\Lu\\L2^n)- (3.8) 

Combining p.7p and p.8p . and dividing through by the weighted iJp-norm of 
u, completes the proof for the case k S [1,2]. 

Now let < K < 1. Again subdivide the domain into i7+ and fl^ by intro- 
ducing a cut along the curve K — 0. Integrate by parts, choosing a = — 1; 

-NIC in n+ 
NIC in n- ' 

where iV is a constant satisfying 

l + i5 l~5 , , 
< < 3.9 

3 - K K + 1 ^ ' 

for a sufficiently small positive constant 5; c = —ANy. The boundary integrals 
involving a and c on either side of the cut will cancel and the boundary integrals 
involving b will be zero along the cut. Inequality (|3.7p can be derived by an ar- 
gument broadly analogous to the case k G [1,2]. Applying (13. Sp then completes 
the proof. □ 

Theorem 3.1 ( [28]). The Dmchlet problem IEM>, P^ . with k £ [0,2] 

possesses a distribution solution u G L^(r2) for every f G H^^{il; IC). 

Proof. Again, we only outline the proof {c.f. [21], Theorem 2.2). Define for 
^ e a linear functional 

Jf{LO^{f,0- 

This functional is bounded on a subspace of by the inequality 

l(/,OI< 11/11 Hl(n-X) (3.10) 

and by applying Lemma 3.1 to the second term on the right. Precisely, J/ is 
a bounded linear functional on the subspace of (fi) consisting of elements 
having the form with ^ € {Vl) . Extending Jj to the closure of this 
subspace by Hahn-Banach arguments, we obtain a functional defined on all 
of . The Riesz Representation Theorem then guarantees the existence of a 
distribution solution in the self-adjoint case. If k 1, then L is not self-adjoint, 
and 

Estimating L for k in [0, 2] will also yield estimates for the adjoint L* . Applying 
the preceding argument to the adjoint operator completes the proof of Theorem 

o □ 
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3.2 Mixed boundary value problems with closed boundary 
data 



It is also possible to form mixed Dirichlet- Neumann problems for operators 
of the form p.4p . Mixed boundary value problems arise in various contexts 
in plasma physics {e.g., [7]) and in related topics from electromagnetic theory 
{e.g., [IH], which is related to the model of Sec. I2.4.I|) . However, the results of 
this section also imply — by taking the set of boundary points on which the 
Dirichlet conditions are imposed to be empty — the existence of weak solutions 
to a class of Neumann problems. 

Denote by u = U2) and w = {wi, 1^2) measurable vector-valued functions 
on Q,. Define 'H.s to be the Hilbert space of measurable functions on 51 for which 
the norm induced in the obvious way by the weighted inner product 



(u,w)^= / / (|/C| 
J Jn 



uiwi + U2W2) dxdy 



is finite. In the notation for these spaces, R denotes a diagonal matrix having 
entries |/C| and 1. 

By a weak solution of a mixed boundary-value problem in this context we 
mean an element u S T-Li\{Vl) such that 

- (u, /:*vif)i2(0;ii2) = (f , w)^2(n.R2) (3.12) 

for every function ^ <E for which Sk-^C*vf E (f2;R2) and for 

which 

wi=Oy{x,y)eG (3.13) 

and 

W2=0y{x,y) edn\G, (3.14) 
where G C dft. Choose the differential operator C to have the form 



ICdx dy \ f K 

dy -d^ J^yo 



(3.15) 



where k is a number in [0, ll 



Theorem 3.2 ( [28 ). Let G be a subset of dH, and let JC — x — . Define the 
functions b {x, y) = mK. + s and c{y) = fiy — t, where fi is a positive constant, 

r {^i + S)/2 tnn+ 

\ {n-S)/2 inQ- 

for a small positive number 5, and t is a positive constant such that fj,y — t < 
OVj/ G ri. Let s be a sufficiently large positive constant. In particular, choose s 
to be so large that the quantities mIC + s, 2cy + s, and b^ + fCc^ are all positive. 
Let 

bdy-cdx<0 (3.16) 
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on G and 

JC {bdy - cdx) > (3.17) 

on d^l\G. Then there exists for every f such that ^^^Ai'^i G L^(r2) a weak 
solution to the mixed boundary-value problem llS.12\) ~ l3.14\ l for C given by Eq. 
i3.15\) with K = 0, where the superscripted T denotes matrix transpose. 

Proof. We give the idea of the proof 28 . One shows that there exists a positive 
constant C such that 

C*M^) J (1^1*? + *2) dxdy 

for any sufficiently smooth 2-vector provided conditions (|3.13l) . (|3.14l) are 
satisfied on the boundary for w = A^*, where C* is given by (I3.15|) with k = 1 
and 

This inequafity leads to an application of the Riesz Representation Theorem by 
arguments which are roughly analogous to those used to prove Theorem l3.1l □ 

Despite the technical nature of the hypotheses in Theorem 13.21 simple do- 
mains which satisfy them are very easy to construct — e.g., a box in the first 
quadrant having a vertex at the origin of coordinates, or a narrow lens about 
the sonic curve in the first quadrant. Note that by taking G to be the empty set, 
we obtain a solution to the closed conormal problem (c./. [32]). But in order for 
Theorem [32] to guarantee a solution to the closed Dirichlet problem, we would 
need to find a domain on which G could be taken to be the entire boundary; 
it is not obvious how to construct such a domain. And, as was also the case 
in Theorem 13.11 the proof of Theorem 13.21 does not establish the uniqueness of 
solutions. 

In addition to its intrinsic mathematical and physical interest, the formula- 
tion of boundary value problems illuminates other topics in the analysis of the 
cold plasma model. For example, it is shown in Sec. 12.4.31 by a tedious ana- 
lytic argument, that away from the origin the governing equation for the model 
is of Tricomi type, whereas in the neighborhood of the origin it is of Keldysh 
type. This distinction is also suggested, without reference to such terminology, 
by other analytic arguments in [S^ and in Sec. 4 of [40]. If we try to form a 
standard elliptic-hyperbolic boundary value problem in which the hyperbolic 
region is composed of intersecting characteristics, we might choose both these 
characteristics to originate at points on the arc of the resonance curve x = y^ 
that lies in the first quadrant, or both of them to lie in the fourth quadrant. 
We then obtain a standard problem for a vertical-ice-cream-cone-shaped region 
(in the former case, the ice-cream cone is held upside down), similar to those 
formulated for the Tricomi equation (Eq. (12.441) with IC{£,,ri) = £,). The domain 
geometry is exactly analogous to, for example. Fig. 2 of [25], with the line AB 
in that figure replaced by an arc of the curve x ^ y'^, lying either completely 
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above or completely below the a;-axis. But the origin will not be included, as 
that is a singular point of the characteristic equation (|2.74|) . If we include the 
origin, we are led to a hyperbolic region bounded by characteristics in the second 
and third quadrants, a horizontal-ice-cream-cone-shaped region similar to those 
formulated for the Cinquini-Cibrario equation (Eq. (I2.46|) with K. (^, rf) = £). In 
this case typical domain geometry is analogous to Fig. 2 of [T^, with the line 
MN in that figure replaced by an arc of the curve x = which is symmetric 
about the x-axis; see also Remark i) following Corollary 11 of [55]. Thus the 
defining analytic character of the equation is clearly apparent in the geometry 
of the natural boundary value problems. 
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